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ABSTRACT 


In this dissertation an attempt has been made to solve a problem of vehicle routing and 
scheduling in a multi-product distributed manufacturing environment. It is being assumed 
that processing facilities are distributed geographically and the material has to be 
transported between these processing centres according to their product processing 
sequence or alternatively transportation between process centres are significant in 
comparison for process time. 

To solve this problem, we formulate the problem as a mixed integer program. As the 
problem size increases it is computationally difficult to solve the problem with exact 
algorithms and the problem is NP-hard. Given the intrinsic difficulty of these class of 
problems, heuristic methods seem to offer the most promise for practical size problems 
We propose two types of heuristic algorithms to solve this problem based on sequential 
and parallel insertion of the tasks. We also test a few despatch rules in assigning the tasks 
to the vehicles. The performance of heuristics is compared computationally. The test 
problem set includes routing and scheduling environments that differ in terms of the type 
of data to generate the problems, the percentage utilisation of process centres, slack time 
available for the products, processing time to travelling time ratio and the problem size 
parameters. Heuristics are also modified and performance studied to situations where the 
number of vehicles is restricted, and where the violation of due-dates of the products are 
allowed by the inclusion of delay costs. 
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CHAPTER 1 

INTRODUCTION 

1.1 Vehicle routing and scheduling : 

Transportation comprises a significant fraction of the economy of most developed nations. For 
example, a National Council of physical Distribution Study [ 1978 ] estimates that 
transportation accounts for 15% of the U.S. gross national product. This economic 
importance has motivated both private companies and academic researchers to vigorously 
pursue the use of operations research and management science to improve the efficiency of 
transportation. [ 3 ] 

Various modes of transportation exist including air, rail, ship and motor vehicle. The research 
on transportation has focused on different issues in each mode. In air, the efficient schedule of 
airline crews [Lavoi, Minous & Odier, 1988] has received primary attention while in rail, use 
of large-scale real time computer control systems has dominated the research agenda ( e.g. 
Harker 1990 ). The literature on both motor vehicles ( trucks, school buses and general 
passenger buses ) and ships has focused on a common problem -- the efficient use of a fleet of 
vehicles that must make number of stops to pick up and/or deliver passengers or products. 
The problem requires one to specify which customers/products should be delivered by each 
vehicle and in what order so as to minimise total cost subject to a variety of constraints such as 
vehicle capacity and delivery time restrictions. [ 3 ] 
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The work in this dissertation is motivated by a vehicle routing and scheduling problem in a 
geographically distributed manufacturing system, for an multi - product or jobbing type of 
organisation. This problem resembles the Pickup - Delivery problem with time windows 
( PDPTW ) and Advance Request Dial-A-Ride with Time Windows ( ADARTW ) with 
additional product sequence constraints ( more details in section 2. 1 ). In this dissertation we 
will concentrate on, minimisation of the total transportation cost involved in transporting the 
products between the processing centres according to their product sequence, without 
violating the time windows. 

1.2 Literature review : 

Laporte etal. [ 1994 ], stated in his survey that “Over the past three decades, operation 
researchers interested in vehicle routing and scheduling have emphasised in the development 
of algorithms for pure vehicle routing problems ( Mingozzi, and Toth 1979, Bodin, Golden, 
Assad, and Ball 1983, Laporte 1992 ) and vehicle routing problems with time windows 
( Desrosiers, Sounds 1985 , Solomon 1987 )”. This section provide a brief review of the 
literature available to solve different vehicle routing problems with time windows. 

In this literature survey, we considered two types of problems. One is the vehicle routing and 
scheduling problem with time windows ( VRSPTW ), which defined as follows. A number of 
vehicles are located at a single depot and must serve a number of geographically dispersed 
customers. Each vehicle has a given capacity. Each customer has a given demand and must be 
served within a specified time window. The objective is to minimise the total cost of travel [1]. 
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The second problem type is the pickup and delivery problem with time windows ( PDPTW ). 
Again, there is a single depot, a number of vehicles with given capacities, and a number of 
customers with given demands. Each customer must now be picked up at his origin durin n a 
specified time window, and delivered to his destination during another specified time window. 
The objective is to minimise total travel cost. [ 1 ] 

The special case in which all customers demands are equal is called the dial-a-ride problem 
( DARP ). It arises in transportation systems for the handicapped and the elderly In these 
situations, the temporal constraints imposed by the customers strongly restrict the total vehicle 
load at any point of time, and the capacity constraints are of secondary importance. The cost 
of route combination is a combination of travel time and customer dissatisfaction. [ 1 3 

Christofides, Mingozzi and Toth [ 1981 ] proposed a dynamic programming algorithm to 
solve the single vehicle problems with time windows. Such a dynamic programming approach 
can be useful for several NP-hard routing problems. However, the cardinality of the state 
space is usually exponential in the problem size. The practical use of dynamic programming in 
this context is restricted to state spaces of at most pseudo polynomial size and relatively small 
problem instances. 

Desrosiers, Dumas and Soumis [ 1986 ] give a similar 2n-stage dynamic programming 
algorithm for a capacitated single-vehicle PDPTW. They proposed number of state elimination 
rules to reduce the computational effort. This algorithm can solve real-life problems with up to 


40 customers. 
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Dumas, Desrosiers and Soumis [ 1991 ] proposed an exact algorithm for a multiple vehicle 
pickup-delivery problem with time windows. This algorithm is based on a Dantzig-Wolfe 
decomposition approach embedded into a branch-and-bound search tree. The master problem 
results in the linear relaxation of a set partitioning type model, while feasible routes or columns 
are generated by a sub problem modelled as a shortest path problem with coupling, 
precedence, time window and capacity constraints. The algorithm has been successful in 
solving two real life problems of sizes 19 and 30 requests, respectively. Very good solutions 
were obtained within the exploration of a few nodes in the branch-and-bound tree, typically 
less than 10. The algorithm is appropriate for the PDPTW when the problem solutions require 
an average up to 5 requests per vehicle, i.e., 10 pickup and delivery stops. However it 
becomes ineffective in the dial-a-ride context due to the actual size of practical problems have 
a large number of requests per vehicle. 

Marius M. Solomon [ 1987 ] had applied Clarke and Wright‘s [ 1964 ] Savings heuristic, Time 
oriented Nearest-Neighbour heuristic, Time oriented sweep heuristic ( original sweep heuristic 
was proposed by Gillett and Miller [ 1974 ] ) to VRSPTW. He also proposed an Insertion type 
of heuristic for this problem. Through his extensive computational study for different type of 
problems, he concluded that the insertion type of heuristics proved to be very successful. 

J.J.Jaw, Odoni, Psaraftis and Wilson [ 1986 ] had proposed an algorithm for an Advance 
Request Dial-A-Ride Problem with time windows ( ADARTW ). In this problem, “advance - 
request” means that all the requests are received well before the time of vehicle dispatching. 
The algorithm uses a sequential insertion procedure to assign customers to vehicles and to 
determine a time schedule of pick-ups and deliveries for each vehicle. 
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1.3 Present work : 

In the literature surveyed no work had been done for a PDPTW with additional sequence 
constraint. Hence as stated earlier the work in this dissertation mainly deals with a PDPTW 
with a sequence constraint case in a multi-product distributed manufacturing system. 

In chapter 2 we define the problem formulate this problem as a conventional Pickup-Delivery 
Problem with Time Windows ( PDPTW ) with additional product sequence constraint. In this 
chapter we also describe the Integer programming formulation for this problem. In chapter 3 
two major heuristic approaches proposed for this problem are described. Further, few 
despatch rules for the selection of the task to be assigned to the vehicle are discussed. We also 
describe modified heuristics for the case, when the number of vehicles is restricted. In chapter 
4 we discuss the design of the problems and study the performance of the heuristics and 
despatch rules for various problems designed. Chapter 5 concludes taking a retrospective view 
on the work done in this dissertation and suggests logical extension to the problem analysed 


for further research. 



6 


CHAPTER 2 

PROBLEM DEFINITION , ANALYSIS AND FORMULATION 

2.1 Problem definition : 

Unlike conventional manufacturing systems, in this thesis it is being assumed that processing 
facilities are distributed geographically and the material has to be transported between these 
processing centres. Such an system implies that, travel times are atleast comparable to the 
processing times and hence transportation schedule and times will have impact on the 
completion time of the job. Further travel costs are significant and will have impact over 
production cost. While such a situation will become relevant in truly global manufacturing 
scenario, it is also meaningful situation in ease of conventional plants when batch sizes are 
small and hence the processing times are small compare to travel time. This also may be 
relevant in a JIT implementation for a multi-product situation, where large number of products 
are required to be produced although in small batches. 

In this thesis, we consider a hypothetical distributed manufacturing system, for an multi- 
product or jobbing type of organisation. In this problem, each batch of product has a pre- 
specified definite processing sequence and due-dates. The processing times at different 
processing centres for each batch of product are also pre-specified. After processing a batch of 
product at a processing centre according to the processing sequence of that product they 
should be transferred to the next processing centre for further processing if any. 
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The objective of the transportation system is to perform all these pickup-delivery requirements 
in time ( i.e., not exceeding the specified due-date ), with minimum total travelling cost . We 

will also consider the case, where the number of vehicles are limited. In such a situation, 
attempt will be to minimise travelling and delay cost. 

Here , the following assumptions are taken into account : 

(1) Vehicle depot is at equal distance from all the processing centres . 

(2) Ail the vehicles have equal capacity . 

(3) All the vehicles have uniform travelling speed . 

(4) Loading and Unloading times are negligible (or) are absorbed in travel times. 

(5) More than one batch of product can be simultaneously processed at a processing centre 
(or) scheduling delays are absorbed in process times. 

(6) Production capacity are sufficiently large and hence there are no schedule delays at the 
process centres (or) they are absorbed in the process times. 

2.2 Problem analysis : 

Let there be ‘n p ’ number of batches of products and l n/ number of processing centres in a 
manufacturing system. A product batch 7 ’ is to be processed in ‘m, ’ number of stages. At 
each stage it is to be processed in a specified processing centre. From the processing time, 
travelling time and product batch sequence data, the time windows of all pickup and delivery 
tasks at all stages can be found for each batch of product as follows : 
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Notations : 

EPT(i J) 

EDT(i J) 

LPT(i,j) 

LDT ( i , j ) 

Dd, 

SI, 

dik 

tl(i , k) 

t2 ( M-ik 7 M-ik+1 ) — > Travelling time from the processing centre at kth stage to processing centre 
at (k+l)th stage for ith batch of product . 

nj — » Total number of pickup nodes for i th batch of product. 

n, + = ( m,~ 1 ) 

n , ' — »Total number of delivery nodes for i !h batch of product. 

n,' =(*1,-1) 

Total number of pickup & delivery nodes for batch i = In, 

n p 

Total number of pickup nodes for all n p batches n = I «j 

i=l 

n p 

Total number of delivery nodes for all n p batches n = Y, n t 

i=l 

Total number of pickup & delivery nodes for all n p batches = In 


— > Earliest pickup time for i* batch of product at ju, stage . 
— » Earliest delivery time for i^ batch of product at j* stage . 
-» Latest pickup time for i* batch of product at jth stage . 

— » Latest delivery time for i* batch of product at j* stage . 
— > Due-date for i* batch of product. 

— > Total slack time available for ith batch of product . 

— > Processing centre for i* batch of product at k* stage. 

— >• Processing time for i* batch of product at ko, stage . 
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2.2.1 Calculation of time-windows for pickup-delivery nodes: 


SI, 


n; 


n; 


= Dd , - ( X tl ( i , k ) + £ t2( pik , [i ±+\ ) ) 

k=l k=l 


j-1 




EPT ( i J ) = I tl ( i , k ) + Z t2( p-ik , (ift+i ) 

k = 1 k = 1 


T T 

(1) (2) 


| (1)- Sum of proc. 
| times of previous 
| stages . 

| (2)- Sum of min. 

| travelling times . 


LPT ( i , j ) = EPT ( i , j ) + SI ; 
EDT(i,j) = EPT(i,j) + tl(i J) 
LDT( i ,j ) = EDT( i , j ) + SI i 


After the calculation of time windows, the present problem resembles the standard pickup - 
delivery problem with time windows, with additional constraint that pickup and deliveries have 
to be made in specified order for each product. 


In this problem, the objective is to minimise the total cost involved in performing all the pickup 
- delivery requirements within the specified time windows . 


The present problem can be better illustrated through the following example. 

2.2.2 Example : 

Let there are two batches of products and each product has three operations. Processing 
sequence, processing times, due dates, vehicle capacity, vehicle speed and distance between 
the processing centres are as follows. 
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Product 1 

Product 2 

Sequence 

MCI -+ MC3 -+ MC2 

MC2 -» MCI -+ MC3 

Proc time 

3 2 1 

2 3 2 

Due-date 

15 

18 


Vehicle capacity = 3 unit loads Vehicle speed = 50 kmph 
Distance between processing centres : 



MCI 

MC2 

MC3 

MCI 

0 

75 

125 

MC2 

75 

0 

100 

MC3 

125 

100 

0 


Now, by using the above data we can calculate the time windows for each pickup- 
delivery task as follows : 

Slack available for 1 st batch of product Sli = 15-((3+2+l)+(2.5+2))=15-10.5=4.5 

Similarly, for the second product Sfe = 7 

Time window for the first p-d task of the first product: 

EPT( 1, 1 ) = 3 LPT( 1, 1 ) = 3+4.5 = 7.5 

LPT( 1, 1 ) = 3+2.5 = 5.5 LDT( 1, 1 ) = 5.5+4.5 = 10 

Similarly, for all other pickup-delivery tasks time windows are calculated. Now, the 


transportation requests are as follows: 
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Pickup 


Delivery 

Product 1 


Productl 

MCI 


MC3 

13.7.5) 


( 5.5. 10) 



Pickup 


Delivery 

Product2 


Product2 

MC2 


MCI 

(2. 9) 


(3.5. 10.5) 


Pickup 

Productl 

MC3 

(7.5. 12) 


Delivery 
Productl 
MC2 
(9.5. 14) 



Pickup 


Delivery 

Product2 


Product2 

MCI 


MC3 

(6.5. 13.5) 


( 9, 16 ) 


Now the objective of the transportation system is to meet all these requests at minimum total 
travelling cost, without violating vehicle capacity, time feasibility and product sequence 
constraints. The present problem can be modelled as a network flow problem. In this network 
pij & dij represents pickup and delivery nodes of 1 th product at j* stage. Here, an arc exists 
between two nodes, only if it satisfies time feasibility and product sequence constraints and 
also all the pickup nodes will have an incoming arc from depot and all delivery nodes will have 
an out going arc to the depot. For the example problem the network will be as shown in figure 
2 . 1 . 



~ '^( ^bepoT ^) 


Fig 2.1 : Network Flow Model 


Directed arc from (or) to the depot 

Undirected arc 

Directed arcs between nodes 
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Now, the objective is to find out the optimal routing of the vehicles from the depot, which 
minimises the total transportation cost involved. The constraints in this problem are: 

• a node should be visited by only one vehicle that too only once. 

• a node should be visited within the specified time window. 

• a vehicle can visit a delivery node only after visiting the corresponding pickup node. 

• the vehicle which visits a pickup node ( say p21 ) should visit the corresponding 

delivervnode also ( i.e., d21 ). 

• a vehicle can visit a higher stage node ( say pl2 ) of a product, only after the completion of 
lower stage tasks, of that product ( i.e., pi 1-dll ). 

• the vehicle load should not be more than its capacity. 

2.3 Mathematical formulation of the problem : 

Let there be n pickup - delivery tasks to be performed including all batches of products . Then 
the total number of pickups and deliveries is equal to 2 n . Then the different sets of nodes in 
the network are represented as follows : 

Set of all nodes in the network N ={ 0 , 1 , 2 , . . . , n-1 , n , n+1 , . . . , 2n , 2n+l } 

Here , node ‘O’ and ‘2n+l’ represents the depot . 

Set of all pickup nodes in the network N p ={ 1,2,3, ....,n} 

Set of all delivery nodes in the network N D = { n+1 , n+2 2n } 




Let there be n p batches of products , then 

Set of all pickup nodes for l^batch of product Ni P = { 1 , 2 , . . . , ni } 

Set of all delivery nodes for l st batch of product Ni D = { n+1 , n+2 , . . . ,n+nj } 

Similarly for all other ( k-1) batches of products we can represent the pickup and delivery 
nodes through systematic node numbering . 

et; — > Earliest pickup (or) delivery time at node £ i’ . 
lti — » Latest pickup (or) delivery time at node ‘i’ . 
q; -> Load to be pickedup or delivered at a node 

Notations gnd_ variables used in the mathemeatical formulation : 


x'i j — > Binary variable Vve { 1,2,- • -,n P } , i e { 1,2,. . ,,n } & j e { 1,2,. . ,,n } 
x V j j = 1 when vehicle V travels from node £ i’ to node £ j’ . 

X V i j = 0 Other wise . 

ah — » Actual pickup (or) delivery time at node £ i’ . 

Ld v ; — » Load of vehicle V at node £ i’ . 

Q — » Capacity of a vehicle . 

Cy — » Cost involved in travelling from processing centre i to processing centre j. 

F — > Fixed cost per vehicle. 
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Formulation : 

n 2n 2n n n 

objective function - Min X X X X v j j - C;j + X Xx V 0 j*F 

v=l 1=1 J=1 v=l J=1 

subjected to — 
n 2n 

X X x v lj =l, Vi6{N p uN D ) <T> 

V=1 j=l 

2n 2n 

X x v ij - X x V j ! = 0 , V i e { N p uN d } , v £ [l,n] <D 

j=i j=i 

n 

X x v 0 j < 1 , V € [l,n] <x> 

J=1 

n 

X X V j(2n +l)< 1 , V 6 [l,n] 

J=l 

2n 2n 

X x v ij - X X V j (n+ i) = 0,VieN p ,ve[l,n] <D> 

j=i j=i 

at, + tl ( M-i , Hn+i) < at ( n +i ) V i e N p 

at( n+ i ) + tl ( Hn+i , M-i+1 ) < at( i+ 1 ) V i e N k p , ke [l,np] <HD 

atj+ ( l-X V ij )*M > at; + 12( n; , (ij ) V i , j £ { N P ^N° } & v e [l,n] *CjT^) 

eti< at ; < lti V i € { N p uN d } >C£) 

Ld v j = Ld v i + qj . x V ij V i , j e{ N p uN° } & v e [l,n] <2D 

LdV 0 V v s [l,n] ►CUT) 

Ld v ; < Q V i e N p 

x v ij e { 0 , 1 } V i , j £ N , ve [ l,n] *Cll) 
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Here, C; j represents the cost involved in travelling from node i to node j and F is the fixed cost 
incurred by using a new vehicle from the depot and M is a very large number. The objective 
function is minimisation of the sum of these costs . 

Constraint 1 ensure that only one vehicle visits a node and that too only once. Constraint 2 
ensure that if a vehicle visits a particular node it should leave that node. Constraint 3 ensure 
that a vehicle should leave the depot only once and constraint 4 ensure that a vehicle should 
arrive the depot only once. Constraint 5 ensure that a vehicle which visites the pickup node, 
that should visit the corresponding delivery node also. Constraint 6 ensure that the vehicle 
should load first at the pickup node and then unload at the corresponding delivery node 
Constraint 7 ensure the product sequence constraint ( i e., only after completion of the 
previous pickup-delivery task, the next p-d task for that batch of product can be started ) 
Constraint 8 calculates the actual pickup and delivery times after visiting the node . Constraint 
9 ensure the actual time of visiting is within the time- windows . Constraints JO , 11 & 12 
ensure the capacity constraint of the vehicles . Consraint 13 ensure non-negativity of variables 

In the above integer programming formulation, the estimated number of variables are of 
0(n 3 ), where ‘n’ is the total number of pickup & delivery nodes. The exact details about the 
number of variables are as follows : 

• Number of binary variables ( i.e., X V y ‘s ) = ( fr + n 2 ) / 2 

• Number of vehicle load variables (i.e., Ld v 1 ‘s ) = n 2 /2 

• Number of time variables ( i.e., at; ‘s ) = n 

• Total number of variables = ( n72 + n 2 + n ) 

The total number of constraints in this formulation are ( 8n 2 + 9n ). 
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Marius M. Solomon [ 1987 ] stated that problems with time windows are frome the 
computational complexity perspective quite difficult. Since the VRP is NP hard, by restriction, 
the VRSPTW is NP hard. Further more, Savelsberge [ 1984 ] has shown that even finding a 
feasible solution to the VRSPTW when the number of vehicles is fixed is itself a NP-complete 
problem. Therefore, the development of heuristic algorithms for this problem class are primary 
interest. 

Yvan Dumas et al. [ 1991 ] had proposed an exact algorithm to solve this type of pickup - 
delivery problem with only time window constraint ( PDPTW ). In their algorithm, they used a 
column generation scheme, with a constrained shortest path as a sub problem. Their algorithm 
worked well for tight capacity and loose time window constraints. The maximum size of the 
problem they solved was having about 50 transportation requests. 

In spite of recent success of optimisation algorithms for vehicle routing with time windows, it 
is unlikely that they will be able to solve large scale practical problems. In many situations one 
has to settle for algorithms that run fast but may produce suboptimal solutions. So, in this 
dissertation we made an attempt to develop heuristics for this problem, which can give 
suboptimal solutions for large scale problems. 
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CHAPTER 3 

SEQUENTIAL AND PARALLEL INSERTION HEURISTIC 

ALGORITHMS 

In this chapter we will discuss two heuristic methods to solve the problem outlined in the 
previous chapter. The two algorithms differ from each other primarily in the fashion in which 
tasks are assigned to the vehicles. In the first approach, the tasks are continued to be assigned 
to a vehicle till no task can be further assigned to it. Only after that the next vehicle is picked 
up for loading. Sequential insertion algorithm uses this approach. In the second approach, the 
tasks are assigned to the best suited vehicle among the vehicles available. Parallel insertion 
algorithm uses this approach. In the section 3.1 we will describe in detail the working of 
Sequential Insertion algorithm. In this section we will also describe all the common features 
for both the algorithms. In the section 3.2 we will describe the Parallel Insertion algorithm. In 
the section 3.3 we will describe the heuristics with simple despatch rules applied for both the 
Sequential and Parallel insertion case for the task selection. In the last section 3.4 we will 
describe the algorithms to consider the case where number of vehicles is restricted. 

Before describing the algorithms, we will define some of the terms used. 

Definitions: 

1) Pickup-delivery task ( p-d task ) : 

The activity of loading a batch of product at a process centre and transporting those 
products to another process centre and unloading that batch is called a pickup-delivery 


task. 
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2) Sequential Feasible tasks : 

Those tasks for which all the preceding tasks have been completed are called sequence 
feasible tasks. 

3) Feasible insertion combination of a task : 

The insertion combination of the pickup and delivery nodes of a task is feasible only if 

(a) it is sequence feasible 

(b) it does not exceed vehicle capacity 

(c) it does not exceed the latest delivery time of that task and also the latest 
delivery times of the tasks which are already in the travel schedule of a vehicle 

4) Feasible task : 

A task is said to be feasible task to assign, to a travel schedule of a vehicle only if it 
has at least one feasible task insertion combination. 

3.1 SEQUENTIAL INSERTION ALGORITHM ( SEQIA ): 

3.1.1 Overview of the algorithm : 

In this algorithm pickup-delivery tasks are assigned to the travel schedule of 
the first vehicle till there is no feasible task, which can be assigned to that vehicle. After this 
the algorithm starts assigning tasks to the next vehicle. Thus this algorithm proceeds, by 
sequentially assigning tasks to the vehicles, until all the transportation requests are processed. 
Central to the algorithm are: 

(i) a search for feasible insertions of tasks into travel schedules 

(ii) an optimisation step to find the best feasible insertion 



19 


• Insertion step: 

An insertion of a task into the travel-schedule is feasible only if it does not lead to the 
violation of vehicle capacity constraint , time feasibility constraint and process 
sequence constraint for the newly assigned task and for all other tasks already assigned 
to that vehicle. 

• Optimisation step: 

The optimisation step deals with minimising the incremental cost due to insertion of the 
new task into a vehicle travel-schedule. The cost function used is the sum of tangible 
cost due to increase in travel distance and intangible costs due to the increase in slack 
utilisation and waiting time( more details in section 3.1.4 ) 

Step by step procedure of this algorithm is illustrated through flowchart in figure 3.1.1 
Notations : 

FTS -- Feasible task set. 

It stores all the process sequence feasible tasks details at every stage of the 
algorithm. For each task it stores the earliest and latest pickup-delivery times 
and also the pickup and delivery process centre numbers. 

VAS( k ) — Vehicle activity set for vehicle ‘k\ 

It stores the details of the travel schedule of vehicle ‘k’ at all stages of the 
algorithm. For each task the following details are stored: 

> Actual time of occurrence of that task. 

> Process centre at which that task occurs. 

> Maximum possible shift 1 for that task. 
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Pit 

Dfc 

Prt ( i, k ) 
Trt ( di,d 2 ) 

Dist ( di, d 2 ) 
Ctr 

Csl 

Cwt 

Cpen 

FC 

Slu( i ) 

PSlu ( 1 ) & E>Slu( i ) 


Sla ( i j ) 

EPT( i j ) 

EDT( i j ) 

LPT(ij) 

LDT(ij) 

APT(ij) 

ADT(ij) 

Pshift 

Dshift 


" Total number of the products. 

-- Pickup process centre number of k* stage task of i* product. 

-- Delivery process centre number of k* stage task of i* product. 

” Processing time of k* stage task of i* product. 

— Travelling time from process centre di to process centre d 2 . 

— Distance from process centre di to process centre d 2 . 

— Travelling cost per unit distance travel. 

— Slack utilisation cost per unit slack time utilised for a product. 

— Waiting cost per unit waiting time utilisation. 

— Penalty cost per unit time violation of due-date for a product 

— Fixed cost per vehicle. 

— Product i’s total slack time utilised due to an insertion of new task 

into the travel schedule. 

-- Product i’s slack time utilised due to new pickup node and delivery 
node insertion respectively 
-- Slack available for product T at jth stage 
-- Earliest pickup time of product i at jth stage 
-- Earliest delivery time of product i at jth stage 

— Latest pickup time of product i at jth stage 

— Latest pickup time of product i at jth stage 

— Actual pickup time of product i at jth stage 

— Actual delivery time of product i at jth stage 

~ Shift due to the insertion of the pickup node ( Local variable ). 

— Shift due to the insertion of the delivery node( Local variable ). 
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Fig 3.1.1 : Flow chart for Sequential Insertion Algorithm 














Fig 3.1.1 : Flow chart for Sequential Insertion Algorithm_( Continue ) 








The sequential insertion algorithm (SEQIA) first calculates the time windows of all the 
sequence feasible tasks ( Refer section 3.1.2 ). Then it stores the details of these tasks in FTS. 
SEQIA start by indexing the task having minimum EPT ( Earliest Pickup Time ) in FTS to the 
first vehicle. Then it determines the actual times of the inserted task.( Refer section 3.1.3 ). 
Then updates the FTS. i.e., it deletes details of the task assigned to the vehicle and add 
the details of the succeeding task in the corresponding product sequence ( Refer section 
3.1.4) Then it processes each task in FTS. Here, each task has a pickup and delivery 
node 1 Now, it starts checking the feasibility of insertion of each task in all possible 
insertion positions in the present travel schedule in a systematic manner (more details in 
Section 3.1.5). Among all the feasible insertion combinations it stores the details of the 
best combination ( i.e., having minimum total incremental cost ) for each task. Then it 
selects the task combination, which has the least incremental cost among all the best 
combinations stored in the previous step( Refer section 3.1.6 ). This task combination is 
actually assigned to the corresponding vehicle travel schedule. The algorithm calculates the 
actual time of occurrence of the newly assigned pickup and delivery tasks and the maximum 
shift possible for those tasks. It also modify the actual time and maximum shift possible values 
for all the nodes which are effected in the travel schedule of that vehicle ( more details in 
Section 3.1.2. ) due to the newly assigned task. It also updates the FTS ( more details in 
section 3.1.3 ). SEQIA repeats this process until all the feasible task insertions are 
completed for that vehicle. If there is no feasible task to insert to that vehicle, then it select 
next vehicle to schedule and repeat the procedure till all the tasks to be done are completed ( 
i.e., till FTS becomes empty ). 


1 Here, a node represents either a pickup or delivery tasks in the schedule . 



3.1.2. Determination of FTS : 

This section describe the detailed calculations involved in the determination of FTS. At the 
beginning, the algorithm determines the time windows for each of the first stage tasks of all 
the products as follows: 

EPT(i, l) = Prt(i, 1); 

LPT ( i, 1 ) = EPT( i, 1 ) + Sla(i, 1 ); 

EDT(i,l) = EPT(i,l) + Trt(P ll ,D ll ); 

LDT ( i, 1 ) = EDT ( i, 1 ) + Sla ( i, 1 ) ; 

3.1.3 Determination and modification of actual pickup and delivery times: 


This section describes the calculation and modification procedure of actual pickup and delivery 
times of tasks. The calculation procedure can be better explained through an example. 
Consider a case as shown in figure 3.1.2. 
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Already scheduled tasks on a vehicle 


Fig 3.1.2:Hlustration of determination and modification of actual pickup and deliveiy 

times 


Here, P^ D« represent the pickup and delivery ( p-d ) nodes of I th product batch at X th stage . 
Similarly, Pj y , Dj y , Pkz & D^z also represent the p-d nodes. Figure 3.1.2 shows that X th and y* 


Stage tasks of i* and j* products respectively, are already scheduled on a vehicle. A new z* 
stage task of k* product is decided to be inserted in the combination as shown in figure 3.1.2. 
The algorithm first calculates the actual pickup time and the shift due to the insertion of the 
pickup node and then it will calculate the actual delivery time and modify the actual 
completion times of the tasks following these nodes by calculating the shift due to the insertion 
of delivery node as follows: 

• Calculation of the actual pickup and delivery times of the newly inserted task and shifts 
due to the insertion : 

APT ( k, z ) = max. { EPT( k, z ) , (APT ( i, x ) + Trt ( Pjx , Pkz ) ) } 


Pshift = max. {0 , ( EPT( k, z) - (APT( i, x) + Trt (P» , P fe )))}+ Trt ( P fa , P jy )-Trt( P ft , P ta ) 


Waiting time 


ADT ( k, z ) = ADT ( i, x ) + Pshift + Trt ( D* , ) 

Dshift = Trt ( Dix , D b ) + Trt ( Db , D jy ) - Trt ( Dix , D jy ) 

• Modification of actual times of the tasks which were effected by the newly inserted task : 

Due to the insertion of the new task, the tasks which are in between the newly 
assigned pickup node and delivery node are shifted by the amount Pshift. The tasks following 
the delivery node get shifted by the amount ( Pshift+Dshift ). 

In the example shown in figure 3.1.2 APT( j, y ) and ADT ( i, x ) get shifted by Pshift 
and ADT ( j, y ) get shifted by ( Pshift + Dshift ). 




3.1.4 Modification of FTS : 


This section describe the detailed calculations involved in the modification of FTS. Whenever 
a task is assigned to a travel schedule, the data in FTS is to be modified accordingly as 
follows: 

• Calculates the time windows for the task succeeding the assigned task in its process 
sequence, as follows: 


EPT ( i, k ) = ADT ( i, k-1 )■ + Trt ( , P* ) ; 

LPT ( i, k ) = EPT ( i, k ) ■ + Sla ( i, k ) ; 

EDT(i,k) = EPT(i,k) + Trt(P ik ,D ik ); 

LDT ( i, k ) = EDT ( i, k ) + Sla ( i, k ) ; 

• Stores the pickup and delivery process centre numbers of the succeeding task. 

• Deletes the information about the assigned task and add the information of the succeeding 
task. 

Due to the newly inserted task, some of the tasks already assigned to the vehicle may be 
shifted. So, the algorithm modifies the EPT ‘s and EDT ‘s of all the effected tasks in FTS. The 
calculation procedure is better illustrated through the example discussed in the previous 
section ( refer figure 3.1.2 ). As explained in the previous section, the algorithm first calculates 
the shift due to the pickup node insertion ( i.e., Pshift ) at the specified position. Due to this 
shift, completion times of all the succeeding nodes are shifted by this amount (i.e., for Pj y , Dix 

& Djy). 



This ‘Pshift is added to the EPT ‘s and EDT ‘s of the effected product tasks ( here i* and j* 
product tasks ) in FTS. Then similarly, the algorithm calculates the shift due to delivery node 
insertion. This ‘Dshift’ is also added to the EPT ‘s and EDT ‘s of the effected product tasks 
( here j* product task ) in FTS. 

3.1.5 Feasibility checking and virtual insertion procedure : 


In the section 3.1.1, it was mentioned that the algorithm finds out the best feasible 
combination for each task in FTS, by inserting p-d nodes virtually in all feasible insertion 
positions. This section will describe the detailed procedure, for checking all the feasible 
combinations for a task selected from FTS in the travel schedule of a vehicle. 

.After selecting a task from FTS, the algorithm starts searching from the last node of the travel 
schedule of the vehicle, for the delivery node of the preceding task, in the process sequence 
of the selected task. If the preceding task was not assigned to that vehicle, or if the task to be 
assigned is not having any preceding task (i.e., it is the first task of the product ), then the 
algorithm starts insertion after the first pickup node of the travel schedule. Otherwise, it starts 
insertion after the delivery node of that preceding task so that process sequence requirement is 
satisfied. 

The insertion procedure can be better illustrated through an example. Consider a travel 
schedule shown in Figure. 3.1.3. Let' the new task to be inserted is ( Pix+i , D^ +1 ). Three 
product tasks at different stages are already scheduled on the vehicle. As described above the 
algorithm first searches for the first insertion position and it starts insertion after the Dix node 
as shown in figure.3.1.3. to satisfy the process sequence constraint. Before inserting the 
pickup node it will check the following : 




First insertion 
position for the 
pickup node 






j 


Already scheduled tasks on a vehicle 

Fig 3.1.3: Illustration of the first insertion position for the pickup node 


• First it will check the vehicle capacity constraint, by virtually adding increase in load due to 
the p-node insertion at the specified position. If the total load at this position does not 
exceed the vehicle capacity, then it further checks the time feasibility constraint, as 
described in the next point. Otherwise the insertion position is moved to the next position 
( i.e., after the node Dj y ). 

• To check the time feasibility constraint, the algorithm first find out the minimum shift 
permissible, by finding the least shift allowed among all the nodes following the 
insertion position. Then it calculates the actual shift to be done for that insertion. If this 
shift is less that or equal to the permissible shift, then the pickup node insertion at that 
position is feasible. Otherwise, the insertion position is moved to the next one. 

After deciding the first insertion position of the pickup node, the algorithm starts checking the 
insertion of delivery node after the pickup node as shown in the figure.3.1.4. Similar to the 
procedure explained for pickup node insertion, the algorithm checks the time feasibility of the 
delivery node insertion, by comparing the permissible shift of the nodes following the insertion 
position of the delivery node. If it is feasible, the algorithm calculates and stores the total 
incremental cost ( Section 3.1.6 ) involved in this combination and also stores the details of 


the combination. 
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Already scheduled tasks on a vehicle 

Fig 3.1.4 Illustration of the first insertion position for the delivery node 


Then move the delivery task to the next insertion and repeat the procedure. If it is not feasible, 
it will simply move the delivery node to the next insertion position and repeat the procedure 
Like this the algorithm stores the details of all the feasible combinations, till the delivery task 
reaches the end of the schedule as shown in the figure. 3. 1.5. 
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Fig 3.1.5 : Illustration of last insertion position for the delivery node 


After finding all the feasible insertion combinations by keeping pickup node at first insertion 
position, the algorithm moves the pickup node to the next insertion position and finds out all 
the feasible insertion combinations as explained above. Repeat this procedure till the pickup 
node reaches the last insertion place. The last insertion combination is shown in the figure 


3.1.6. 




Pix+i ‘ \ Dix+i ; 



Already scheduled tasks on a vehicle 
Fig 3.1.6 : Illustration of the last insertion combination of a p-d task 


After checking all the combinations of that task , the algorithm selects the combination having 
minimum total incremental cost. It retains the details of that combination and the cost involved 
in that. Then it deletes the other combination details. After this the algorithm selects the next 
task in FTS and repeat the same procedure. 

3.1.6 The optimisation procedure : 

Inorder to select the best feasible insertion combination, among all the feasible insertion 
combinations of a task into the travel-schedules of the available vehicles, we use an objective 
function designed to evaluate the incremental cost involved in each insertion. 

This cost is taken to be, the sum of the tangible and intangible costs involved in an insertion. 
Here, operating cost involved in travelling is considered as tangible cost. On the other hand to 
control the utilisation of the slack available for each product, slack utilisation cost and vehicle 
waiting time cost are considered as intangible costs. 

The algorithm computes the total incremental cost for each feasible insertion combination. 
This computation procedure can be illustrated through the following example. Consider a 
feasible insertion combination as shown in figure3.1.7. in a vehicle travel schedule. 
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Fig 3.1.7 Dlustration of the incremental cost calculation 

The total incremental cost involved in the insertion is calculated as follows. 

.p . . . roS t + slack utilisation cost 

Total cost = Total increase in travel cost + vehicle waiting cos'- 

Total increase in travel cost = Increase in travel cost due to pickup node insertion 

, . a t n delivery node insertion 

+ Increase in travel cost due to j 

Increase in travel cost due to p-node insertion is given by 

Ctr . [ Dist (P* , Pb) + Dist (Pb , Pjy) ' DlSt (P “ ’ Pjy) ] 

Increase in travel cost due to d-node insertion is given by 

Ctr » [ Dist (Dix , Dt) + Dist (D ta , Pjy) ' DlSt ^ ’ ° jy) 1 

Vehicle waiting cost = Cwt * [ max. { 0 , [ EPTkz - ( APT tf + Trt( tx, kz))]}] 

Up 

Slack utilisation cost = Csl * [ 2 Slu ( i ) ] 

i=l 

, , ... A ( slu ) in a feasible insertion into 

Here, the algorithm calculates all the products slack utilised t 

. . Th ; s calculation procedure can be 

the travel schedule to calculate the slack utilisation cost. 

better illustrated through the example shown in figure 3.1-7- 

, «. oickup node is virtually inserted 

It first calculates the slack utilised for the product, whose f 

r_ii 0 wing the pickup node insertion 
(i.e., here, k* product ) and the products whose tasks are i° u 

position (i.e., here I th and j* products ) as follows: 



Product k* slack utilised PSlu ( k ) = APT k 2 - EPT kz ; 

Products i and j’s slack utilised PSlu ( i ) & PSlu ( j ) 

= max. { 0,( EPT kz - ( APT ix + Trt (P^)} + Trt (P„ , Pkz) + Trt (P te , P jy ) - Trt (P» , P jy ) ; 
For all other products PSlu is assigned to be zero. 

Then the slack utilised in the delivery node insertion is calculated, similarly as described above 
for pickup node insertion as follows: 

Product k* slack utilised DSlu ( k ) = (ADT kz - EDT kz ) - PSlu ( k ); 

Product j* slack utilised DSlu ( k ) = Trt (D* , D te ) + Trt (D^ , D jy ) - Trt (D* , D jy ) ; 

For all other products the DSlu is assigned to be zero. 

The total slack time utilised for all products ( Slu ) in the insertion is sum of PSlu and DSlu. 
Here, the performance of the algorithm depends upon the values of C st and Cwt . According to 
these values the slack utilisation is controlled. The optimum values of these parameters are 
problem dependent. So, here the algorithm is repeated for various values of C st and Cwt 
, within a range. It stores the details of the schedules having minimum total cost, including the 
fixed cost associated per vehicle. 


3.2 PARALLEL INSERTION ALGORITHM ( PARIA ): 

3.2.1 Overview of the algorithm : 

The major difference between PARIA and SEQIA is that, in PARIA inserting 
on a new vehicle is started whenever a task has no feasible insertion combination on all the 



existing % e hides. Where as in SEQIA, insertion on a new vehicle is started, only when there 
are no feasible tasks to insert on all the existing vehicles. In PARIA a task is inserted on all the 
available vehicle travel schedules, in all feasible combinations and select the best one among 
those is selected. Similar to the SEQIA, an insertion combination of a task into a travel 
schedule is feasible only if it does not lead to the violation of vehicle capacity, time feasibility 
and process sequence constraints, for the newly assigned task and for all other tasks already 
scheduled on the vehicles. The optimisation step deals with minimising the incremental cost 
due to insertion of new task into a vehicle travel schedule. 

Step by step procedure of the algorithm is illustrated through flow chart in figure 3.2.1. 

The parallel insertion algorithm ( PARIA ), first calculates the time windows to all the 
sequence feasible tasks. Then it stores the details of these tasks in FTS. Similar to SEQIA, this 
algorithm also starts by indexing the task having minimum EPT in FTS to the first vehicle. 
Then update the FTS and starts insertion procedure for the new task assignment on the first 
vehicle. This process of insertion on the first vehicle is same as in SEQIA, till for a sequence 
feasible task no feasible insertion combination is on this vehicle i.e., task is infeasible on this 
vehicle. Then the algorithm insert this first such task on the second vehicle. Now, the 
algorithm starts checking for the feasible insertion combinations of all the tasks in FTS, on 
both the vehicle travel schedules. It selects the task combination which has minimum 
incremental cost, among all the feasible insertion combinations. The virtual insertion procedure 
in PARIA is similar to SEQIA, except that it has to check for the best feasible insertion on all 
the available vehicle travel schedules. The optimisation procedure is also same as in SEQIA. 
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Fig 3.2.1 : Flow chart for Parallel Insertioa Algorithm 












© 




^ 1 171 


,1 r — TncAH-mn Alanrithm ( Continue ) 












3.3 HEURISTICS WITH DESPATCH RULES : 


In both the heuristic approaches, the task selection is done by optimising the incremental cost 
with the various values of slack cost and vehicle waiting cost. This approach is quite time 
consuming and computationally intensive. In this section we describe four heuristics based on 
simple despatching rules for task selection. These rules are used for both sequential and 
parallel insertion heuristics. 

3.3.1 Earliest due date procedure ( EDD ) : 

Jackson [ 1955 ] proved that scheduling according to EDD rule, will minimise the maximum 
lateness. As we are looking for a travel schedule in which all tasks are to be completed within 
their due-dates, selection of the task for insertion using EDD is likely to be effective. 
According to this rule, the algorithm selects the task which has the minimum latest delivery 
time. After selecting the task, the algorithm virtually inserts the pickup and delivery nodes of 
this task in all feasible insertion positions. If there is no feasible insertion for this task, then it 
selects the next best according to this rule. Among all the virtual insertion combinations of the 
selected task, the algorithm selects the combination which has the least incremental travelling 
cost. 

3.3.2 First come first serve ( FCFS ) : 

According to this rule the algorithm selects the task which has minimum earliest pickup time 
and virtually inserts the pickup and delivery nodes of this task in all the feasible insertion 
positions. If there is no feasible insertion for this task, then it selects the next best task 



according to FCFS rule. Among all the virtual insertion combinations of the selected task it 
chooses the combination which has the minimum incremental travelling cost. The motivation 
behind this rule is that, the vehicle is loaded with a task as soon as the task is available. This is 
likely to result in reduced waiting times for the vehicles. 

3.3.3 Minimum slack per operation ( MINSL ): 

According to this rule, the algorithm selects the task of a product which has minimum slack 
available per operation at that stage. After the selection, it virtually inserts the pickup and 
delivery nodes of this task in all the feasible insertion positions. Among all the virtual insertion 
combinations of this task, it chooses the combination which has the least incremental cost 
This rule will effectively give priority to those tasks which have smaller slack available. 

3.3.4 Systematic slack costing ( SYSSL ) : 

In this heuristic penalty is imposed for over utilisation of the slack per operation available for a 
task. For a task insertion combination, difference between actual slack utilised and the slack 
available per operation for that product are computed. Slack over utilisation cost is computed 
as the product of penalty per unit slack over utilisation and the excess slack utilised. We have 
selected the penalty for over utilisation of the slack equal to the travel cost per period. 

3.4 Restricted number of vehicles at the depot : 

The algorithms discussed in the previous sections deal with the problems for which the 
violation of due-dates ( i.e., hard time windows ) are not allowed and also without any 
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restriction on the number of vehicles available at the depot. But, in the practical scenario, the 
number of vehicles available at the depot are restricted. So, we made an attempt to observe 
the performance of these algorithms in case of restricted number of vehicles at the depot In 
this case the violation of due-dates are allowed, when the available number of vehicles are not 
sufficient to meet all the transportation requests. 

To control the violation of the due-dates of the products, a penalty cost per period violation of 
due-date is included in the calculation of incremental cost in the algorithms. Whenever, a new 
task insertion violates the due-date of that task or the due-dates of the tasks which are already 
scheduled to that vehicle: the algorithm calculates the penalty cost involved in that insertion 
and add this to the incremental cost( Refer section 3.4.1). 

3.4.1 Calculation of penalty cost : 

The calculation procedure can be better illustrated through an example. Consider a sequence 
feasible insertion combination as shown in figure 3.4.1. 



Already scheduled tasks on a vehicle 

Fig 3.4.1: Illustration of calculation procedure of penalty cost 


If this task insertion combination violates the due-dates. The algorithm calculates the penalty 


cost as follows: 
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• if ADT ( k, z ) > LDT ( k, z ) 

then 

tpen ( k, z ) = ADT ( k, z ) - LDT ( k, z ); 

• for the tasks which are between the pickup and delivery node insertion positions ( i.e., Pjy 
and Dix), the algorithm compares the shift due to the pickup node insertion and the shift 
available for those tasks. If the available shift is less than the actual shift due to the pickup 
node ( Pshift ), then the algorithm calculates the difference between these and add this to 
the tpen( k, z ). 

• for the tasks which are following the delivery node, the algorithm compares the shift due 
to both the pickup and delivery nodes insertion and the shift available for those tasks. If 
the available shift is less than the actual shift due to both the pickup and delivery node 
( Pshift+Dshift ) insertion, then the algorithm calculates the difference between these and 
add this to tpen( k, z ). 

Penalty cost = Cpen*tpen( k, z ); 
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CHAPTER 4 

COMPUTATIONAL STUDY 

In this chapter we will describe the details of computational study conducted for both the 
sequential and Parallel Insertion algorithms and the dispatch rules. In the section 4.1 we will 
describe the design of the test problems, which are used for the performance analysis of the 
heuristics. In section 4.2 we will describe the performance analysis of all the algorithms 
discussed in the previous chapter. In section 4.3 we will study time complexity of algorithms 
In section 4.4 we will discuss the performance of heuristics, when number of vehicles is 
fixed. 

4.1 Design of test problems : 

A variety of test problems were generated randomly. In this section we will describe the 
details of the generation of these test problems. In the design of test problems several factors 
that can effect the behaviour of heuristics were considered. We divided these factors into 
following two categories; 

I . Problem process data : 

( 1 ) Processing centres utilisation, 

(2) Slack available for the products, 

( 3 ) Processing time to travelling time ratio. 

II . Problem size : 


( 1 ) Number of products, 

(2) Number of operations per product, 
( 3 ) Number of processing centres. 
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To generate the problems of required class, we programmed a problem generator. The 
complete details of the working of this generator are illustrated through flow chart as shown 
in figure 4.1. In this, we used a standard routine available to generate the pseudo random 
number. The first seed to this random number generator is given through standard input. 

4.1.1 Process data variation: 

In this section we will describe the procedure followed to generate the problems by varying 
the process data of the problem. We generated twenty seven problems, by keeping the 
factors- process centre utilisation, slack time available and processing time to travel time 
ratio at three levels. The actual combinations used are described in Appendix-A. The details 
of these factor levels are discussed below: 
i) Process centre utilisation : 

For this factor the following three levels are considered. 

Level -1 : Equal utilisation of all the processing centres. 

For example if there are four process centres. Each process centre has 
25% utilisation. 

Level -2 : A combination of high, medium and low process centres 
utilisation. 

For example if there are six process centres. Two of them will have 
high percentage utilisation ( say 35%, 30% ), the other two will have 
medium percentage utilisation ( say 10%, 15% ) and the last two 


problems will have low level of utilisation ( say 3%, 7% ). 
Level -3 : A combination of high and low process centres utilisation. 



START 


INPUT 


Average number of products in the problem.( N ) 
Average number of operations per product. ( 0 ) 
Average number of process centres. ( P ) 
Processing time to travelling time ratio. ( kl ) 
Due-date to total processing time ratio. ( k 2 ) 


Generate a random number using uniform distribution in the 
interval [ 0.75N,1.25N ] and consider this as actual number of 
products ( n ) in the problem. 


Generate a random number using uniform distribution in the 
interval [ 0.75P,1.25P ] and consider this as actual number of 
process centres ( p ) in the problem. 


Input the percentage utilisation for each of the processing centres. 

The cumulative percentage utilisation of each of the process 
centres are calculated . The total sum of %utilisation of all process 
centres is equal to equal to 100. 


Generate for each product a random number using uniform 
distribution in the interval [ 0.75*0, 1.25*0 ],. and assign it to the 
number of operations of that product. 


Fig 4.1 Flow chart for the problem generator 
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Fig 4.1 Flow chart for the problem generator( Continue ) 
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For example if there are four process centres.Two of them have high 
percentage utilisation ( say 42%, 43%), the other two have low level 
utilisation ( say 8%, 7%). 
iil Slack time available per product : 

For this factor the following three levels are considered. 

Level- 1 : High 

Slack time available for a product is equal to the total process time of 
that product. 

Level-2 : Medium 

Slack time available for a product is equal to 50% of the total 
processing time of that product. 

Level-3 : Low 

Slack time available for a product is equal to the 10% of the total 
processing time of that product. 

iii) Processing time to travelling time ratio : 

For this factor the following three levels are considered. 

Level-1 : High 

At this level ratio is kept at 4. It means that comparatively the 
travelling time involved in the problem is 25% of the processing time. 
Level-2 : Medium 

At this level ratio is kept at 2. It means that comparatively the 
travelling time involved in the problem is 50% of the processing time. 
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Level-3 : Low 

At this level ratio is kept at 1. It means that comparatively the 
travelling time involved in the problem is 100% of the processing 
time. 

4.1.2 Problem size variation: 

In this section we will describe the procedure followed to generate the problems by varying 
the size of the problems. We generated forty problems, by keeping the factors— number of 
products, number of operations per product and number of processing centres at two levels 
The actual combinations used are described in Appendix-B. For each combination five 
problems are generated with different process data The details of these factor levels are 
discussed below: 

i) Number of products : 

For this factor the following two levels are considered. 

Level- 1 : High 

At this level number of the products for a problem is generated 
randomly in the range 10-18. 

Level-2 : Low 

At this level number of the products for a problem is generated 
randomly in the range 4-10. 


ii) Number of operations per product : 

For this factor the following two levels are considered. 
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Level- 1 : High 

At this level number of the operations per product for a problem is 
generated randomly in the range 10-20. 

Level-2 : Low 

At this level number of the operations per product for a problem is 
generated randomly in the range 5-10. 
in) Number of process centres : 

For this factor the following two levels are considered. 

Level- 1 : High 

At this level number of process centres for a problem is generated 
randomly in the range 10-15. 

Level-2 : Low 

At this level number of the process centres for a problem is generated 
randomly in the range 5-10. 

The detailed computational study is explained in the next section. 

4.2 Performance analysis : 

In this section we will describe the computational performance of the sequential and parallel 
heuristics discussed in the previous chapter and the performance of the heuristics with 
despatch rules, which were discussed in section 3.3. We programmed all these algorithms in 
CJanguage and performed computational test on all the problem classes described in the 
previous section. All the computational experiments are done on HP9000 platform using 


ITNTX nnprfltina cvc+pm 
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To compare the performance of the algorithms and des patch rales (hem ^ we 

considered the following two measures : 

1) Average % deviation of the total travelling cost of an ~ . 

' b n algorithm from the best solution 

among all the heuristics. 

2) Average number of vehicles used. 

Graphs for different problem sets were plotted by taking the average percentage deviation of 
total travelling cost on the y-axis and the despatch role used on the x-axis (Fig 4.2 & Fig 4 4) 
and also by taking the average number of vehicles used on y.^ and , he despatch ^ ^ 

on the x-axis ( Fig 4.3 & Fig 4.5 ). The details of the behaviour of these algorithms for 
different problem sets are discussed below. 

4.2.1 Processing centre utilisation : 


By keeping the processing centre utilisation at three diSe rem combinations, we performed 
computational test for nine problems corresponding to each of the combination of slack 
available and processing time to travel time ratio. Following i s observed 
( 1 ) Average percentage deviation of total travelling cost^Q^ ^ est Qne 
According to this measure: 

a) at the best values of slack cost and vehicle waiting cost ^ Paral|el 

algorithm is performing better than the Sequential Insertion for ^ ^ ^ 

combinations of processing centres utilisation^ Refer Fig 4 2 1) 

b) Among the despatch rules applied to Sequential h sertion ^ Systematic 

Slack Costing rule is performing better than other three rules in case of the even and 
also Kgh-Medium-Low utilisation combination of the procession centres 



the case of high-low utilisation combination of the processing centres 'First Come 
First Serve is performing better than others. .( Refer Fig 4.2.1) 

c) Among the despatch rules applied to Parallel Insertion algorithm the first come 
first serve rule is performing better than the other rules for all the three combinations. 

d) Among all the despatch rules applied for both the algorithms, the first come first 
serve applied to Parallel Insertion procedure is performing better ( Refer Fig 4.2. 1 & 
Table 4.3 ). 

(2) Average number of vehicles used : 

According to this measure: 

(a) Both the sequential and Parallel Insertion algorithms are performing almost 
equally at the best values of slack utilisation cost and vehicle waiting cost. 

(b) For both the algorithms, the First Come First Serve rule is performing better than 
the other rules. 

(c) The First Come First Serve rule applied to the Parallel Insertion case is performing 
better among all the despatch rules applied to both the algorithms for all the three 
combinations. (Refer Fig 4.3) 

4.2.2 Processing time to travelling time ratio : 

By keeping processing time to travelling time ratio at high, medium and low levels, the 
computational tests are performed for nine problems corresponding to each of the 
combination of process centre utilisation and slack available. The details of these results are 


described as follows. 



49 


( 1 ) Average percentage deviation of total travelling, cost : 

According to this measure: 

(a) For all the three levels of processing time to travelling time ratio, the parallel 
insertion algorithm is performing better than the Sequential Insertion algorithm at the 
best values of slack cost and vehicle waiting cost. 

(b) Among the despatch rules applied to Sequential Insertion algorithm the first come 
first serve rule is performing better than the other three rules when this ratio is at high 
level. Systematic Slack Costing is performing better when this ratio is at medium 
level and Earliest Due Date rule is performing better when this ratio is at low level. 

(c) Among all the despatch rules applied to the Parallel Insertion algorithms, the 
First Come First Serve rule is performing better than others at all the three levels. 

(d) Among all the despatch rules applied for both the algorithms, the First Come First 
Serve applied to Parallel Insertion procedure is performing better ( Refer Fig 4.2 & 
Table 4.3 ). 

( 2 ) Average number of vehicles used : 

According to this measure 

(a) Both the sequential and Parallel Insertion algorithms are performing almost 
equally at the best values of slack utilisation cost and vehicle waiting cost. 

(b) For both the algorithms, the First Come First Serve rule is performing better than 
the other rules. 

(c) The First Come First Serve rule applied to the Parallel Insertion case is 
performing better among all the despatch rules applied to both the algorithms at all 
three levels (Refer Fig 4.3 & Table 4.4 ). 
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4.2.3 Slack time available per product : 

By keeping the slack time available per product at high, medium and low levels for all the 
products, the computational test are performed for nine problems corresponding to each of 
the combination of process centre utilisation and processing time to travel time ratio. The 
details of these results are described as follows. 

( 1 ) Average percentage deviation of total travelling, cost : 

According to this measure 

(a) For all the three levels of slack time available per product, the Parallel Insertion 
algorithm is performing better than the Sequential Insertion algorithm at the best 
values of slack cost and vehicle waiting cost. 

(b) Among the despatch rules applied to the Sequential Insertion algorithm the first 
come first serve rule is performing better than the other three rules when all the 
products in the problem have slack availability at high level or at low level, 

Systematic Slack Costing rule is performing better when all the products in the 
problem have slack availability at medium level. 

(c) Among all the despatch rules applied to the Parallel Insertion algorithms, the 
First Come first Serve rule is performing better than others at all the three levels. 

(d) Among all the despatch rules applied for both the algorithms, the First Come First 
Serve applied to Parallel Insertion procedure is performing better ( Refer Fig 4.2 & 
Table 4.3 ). 

f 2 ) A verage number ofvehicles used: 

According to this measure: 



(a) Both the sequential and Parallel Insertion algorithms are performing almost 

equally at the best values of the slack utilisation cost and vehicle waiting cost, when 
the slack available for all products are at high and medium levels. * 

(b) In case the slack availability is at low level for all the products, the parallel 
insertion method is performing better than the sequential insertion. 

(c) For both these algorithms when the slack available for all products are at high and 
medium levels, the despatch rule First Come First Serve is performing better than the 
other three rules. 

(d) In case the slack availability is at low level. Earliest Due Date rule is performing 
better for both the algorithms. 

(e) Among all the despatch rules applied to both the algorithms, the first come first 
serve rule applied to the Parallel Insertion algorithm is performing better, when the 
slack available for all products are at high and medium levels(Refer Fig 4.3 & Table 
4.4 ). Similarly, the Earliest Due Date rule applied to the Parallel Insertion algorithm 
is performing better when the slack availability for all products are low. 

4.2.4 Problem size factors: 

As mentioned already, number of products, number of operations per product and number of 

departments are considered as three size parameters. By keeping each of these parameters at 

two levels, we generated six sets of problems. For each problem set we generated twenty 

problems. The details of these results ar^.a s follows : 

% 

( 1 ) Average percentage deviation of total travelling- cost : - 

CENTRAL* -’jHRAR* 

<-**• A 123288 


According to this measure 
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(a) For all the six set of problems, the Parallel Insertion algorithm is performing 
better than the Sequential Insertion algorithm at the best values of slack cost and 
vehicle waiting cost. 

(b) Among the despatch rules applied to the Sequential Insertion algorithms. 
Systematic Slack Costing rule is performing better than all other rules for all the six 
set of problems. 

(c) Similarly, among all the despatch rules applied to the Parallel Insertion 
algorithms, the First Come First Serve rule is performing better than others except in 
the case of number of products are low, in this case Earliest Due Date despatch rule 
is performing better. 

(d) Among all the despatch rules applied for both the algorithms, the first come first 
serve applied to Parallel Insertion procedure is performing better, when the number 
of products are at high level, number of operations per product are low and number 
of departments are low. The Systematic Slack Costing rule is performing better 
when the number of products are low, number of operations per product are high and 
number of departments are low ( Refer Fig 4.4 & Table 4.5 ). 

( 2 ) Average number ofvehicles used : 

According to this measure 

(a) Both the sequential and Parallel Insertion algorithms are performing almost 
equally for all the six sets of problems. 

(b) Among all the despatch rules applied to both the algorithms, the first come first 
serve rule applied to the Parallel Insertion algorithm is performing better for all six 
type of problem sets (Refer Fig 4.4). - 

Summery of this analysis is shown in table 4. 1 and table 4.2. 
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BEST PERFORMED DESPATCH RULES ON DIFFERENT 

PROBLEMS 


( In terms of average%deviation of total travel cost from the best cost) 


FACTOR 

LEVEL 

Best 

performed 
dispatching 
rule in case 
of sequ. 
insertion 
method 

Best 

performed 
dispatching 
rule in case 
of parallel 
insertion 
method 

Among all 
the 

dispatching 
rules for 
bothe the 
algorithm 
best 

performed 

one 


Even 

SYSSL 

FCFS 

FCFS ( Pari ) 

PROCESS CENTRE 




UTILISATION 

High-Medium-low 

SYSSL 

FCFS 

FCFS ( Pari ) 


High-Low 

FCFS 

FCFS 

FCFS ( Pari ) 


High 

FCFS 

FCFS 

FCFS ( Pari ) 

PROCESSING TIME 




TO TRAVELLING 

Medium 

FCFS 

FCFS 

FCFS ( Pari ) 

TIME RATIO 





Low 

SYSSL 

FCFS 

FCFS ( Pari ) 


High 

FCFS 

FCFS 

FCFS ( Pari ) 

SLACK TIME 





AVAILABILITY 

Medium 

SYSSL 

FCFS 

FCFS ( Pari ) 


Low 

FCFS 

FCFS 

FCFS (Pari) 

NUMBER OF 

High 

SYSSL 

FCFS 

FCFS ( Pari ) 

PRODUCTS 





Low 

SYSSL 

FCFS 

SYSSL (Seq) 

NUMBER OF 

High 

SYSSL 

FCFS 

SYSSL(Seq) 

OPERATIONS PER 





PRODUCT 

Low 

SYSSL 

FCFS 

FCFS(Parl) 

NUMBER OF 

High 

SYSSL 

FCFS 

SYSSL(Seq) 

PROCESSING 





CENTRES 

Low 

SYSSL 

FCFS 

FCFS (Pari) 


Table 4.1 
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BEST PERFORMED DESPATCH RULES ON DIFFERENT 

PROBLEMS 


( In terms of average number of vehicles used) 


FACTOR 

LEVEL 

Best 

performed 
dispatching 
rule in case 
of sequ. 
insertion 
method 

Best 

performed 
dispatching 
rule in case 
of parallel 
insertion 
method 

Among all 
the 

dispatching 
rules for 
bothe the 
algorithm 
best 

performed 

one 


Even 

FCFS 

FCFS 

FCFS ( Pari ) 

PROCESS CENTRE 





UTILISATION 

Fligh-Medium-low 

FCFS 

FCFS 

FCFS ( Pari ) 


High-Low 

FCFS 

FCFS 

FCFS ( Pari ) 


High 

FCFS 

FCFS 

FCFS ( Pari ) 

PROCESSING TIME 





TO TRAVELLING 

Medium 

FCFS 

FCFS 

FCFS ( Pari ) 

TIME RATIO 






Low 

FCFS 

FCFS 

FCFS ( Pari ) 


High 

FCFS 

FCFS 

FCFS (Pari) 

SLACK TIME 





AVAILABILITY 

Medium 

FCFS 

FCFS 

FCFS ( Pari ) 


Low 

FCFS 

FCFS 

FCFS (Pari) 

NUMBER OF 

High 

FCFS 

FCFS 

FCFS ( Pari ) 

PRODUCTS 






Low 

EDD 

FCFS 

FCFS ( Pari) 

NUMBER OF 

High 

FCFS 

FCFS 

FCFS(ParI) 

OPERATIONS PER 





PRODUCT 

Low 

FCFS 

FCFS 

FCFS(ParI) 

NUMBER OF 

High 

FCFS 

FCFS 

FCFS(Parl) 

PROCESSING 





CENTRES 

Low 

FCFS 

FCFS 

FCFS (Pari) 


Table 4.2 
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The following is observed : 

%deviatign of total travelling cost <yf_ a heuristic from the best onej 
According to this measure: 

• At the best values of slack cost and vehicle waiting cost Parallel Insertion algorithm is 
performing better than the Sequential Insertion algorithm, irrespective of the problem 
type. 

• Among the despatch rules applied to the Sequential Insertion algorithm, Systematic Slack 
Costing rule is performing better inmost of the cases( 67 % ). In some of the cases First 
Come First Serve rule( 33 % ) is performing better. 

• Among the despatch rules applied to the Parallel Insertion algorithm. First Come First 
Serve is performing better in all the cases. 

• Among the despatch rules applied to both the Sequential and Parallel Insertion 
algorithms the First Come First Serve rule applied to Parallel Insertion algorithm is 
performing better except when the number of products are low, number of operations per 
product are high, and number of process centres are low. In these three cases Systematic 
Slack Utilisation applied to Sequential Insertion algorithm is performing better. 


Average number of vehicles usedj 
According to this measure 

_• First Come First Serve despatch rule is performing better in both the Sequential and 
Parallel Insertion algorithms. 

_• First Come First Serve despatch rule applied to Parallel Insertion algorithm is performing 
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4.3 Time complexity of algorithms : 

After testing the computational performance of the algorithms, we were interested to check 
the performance of the algorithms in terms of computational time taken for different sizes of 
the problems. The fig 4.6 illustrates the growth rate of computational time taken by these 
algorithms with respect to the number of pickup-delivery tasks in the problem. By observing 
the graph we can conclude that the growth rate of the time taken by the Parallel Insertion 
algorithm is more than that of Sequential Insertion algorithm. The worst case analysis of 
these algorithm shows that the Sequential Insertion algorithm has 0(n 3 ) and for Parallel 
Insertion algorithm 0(n 4 ) where ‘n’ is the total number of tasks in the problem. 

4.4 Performance analysis in case of restricted number of vehicles 
available: 

We performed a computational test for five problems and observed the behaviour of all the 
heuristics, by changing the weight of the penalty cost with respect to the travelling cost. 
Similar study is also comducted by changing the number of vehicles available at the depot. 
The computational results of an example problem are described in table 4.5,4.6,4.7 and 4.8. 
The observations made through these results are as follows : 

• At the best values of slack cost and vehicle waiting cost Parallel Insertion algorithm is 
performing better than Sequential Insertion algorithm. 

• Despatch rules applied to Sequential Insertion procedure are performing effectively when 
compared to the despatch rules applied to the Parallel Insertion procedure. 
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• Minimum slack per operation rule is performing better in case of Sequential Insertion 
method, in controlling the violation of due-dates. Where as the Systematic Slack Costing 
rule is performing better in reducing the travelling cost, but it is not able to control the 
violation of due-dates of the products. 

• Earliest due date despatch rule is performing better in case of Parallel Insertion method in 
controlling the violation of due-dates. Similar to the sequential case, the Systematic Slack 
Costing rule is able to reduce the travelling cost, but it is not able to control the violation 
of due-dates of the products. 

• Results shows that, the behaviour of the algorithms and the task despatch rules, are not 
changing by varying the number of vehicles available at the depot . 

• Results shows that, by increasing the weightage for the penalty cost with respect to the 
travelling cost is able to control the violation of due-dates of the products till some value, 
then the algorithm gets saturated and gave the same results by further increase of the 
weightage. 



Comparision of the performance of all heuristics interms of %dcv. of total travelling cost 

with respect to the best one 
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Avg No. Of Vehls 




































Computational perf ormance of all the heuristics interms of travelling cost with respect to the host 
on the problems generated by varying the process data of the problems 
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High 
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Table 4.3 
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Comparision of the performance of all heuristics in terms of%dev. of total travelling 

cost with respect to the best one 












Comparision of the performance of all the heuristics In terms of avg. numberof 

vehicles used 
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Comparision of the performance of all the heuristics interms of total travelling cost with respect to best one 

on the problems generated bv varying size perameters 
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Time comparision of sequential and parallel insertion algorithms 




Sequential insertion procedure for restricted number of vehicles 




js <o vo 

co co Cj to oo 

O ’S 2 ^ ro 

^ vo cn 



</o — * 0 > to 

O "< *t On r, 

o (N n 

oo r- r- r- 


oo 

00 

00 

to 

r~ 

r- 

Ov 

<N 

1 ^ 

to 

•— < 



_ 2 o 
cc cc cs 

o ^ ^ 

H p « )g 


7^ o> w m «o »n 

O co O on — < rj- © 

H i3 u h On 00 0C 


^ 1 !Q *n 

in -h “ r^ ( 
VO rt CO 



re o - 5 N ® ® s s 

n £ o S r 1 n r< x ^ ^ “ 

r° £ g g oe oo oo o oc oo ts 

H Jo u g hn m *n r- — — — . 

50 OO 00 OC 0C ON On ON 


CO 

2) 

w in cn — < 



^ n — 'to 

o 

CO 

R O ro -0“ 

H 

l-H 

U CN Ov Os 


*oj w H k on 

«g > £ O 

V Tf vo 

o CO O V 

1 co vo 

H £3 ° cx 

5 On On 



£> £ *5 

J2 ^ £ CO 

CO Jx£ O 
> O „£> O 
W CO O 

S3 ^ 2 

m jg 





















Parallel insertion procedure for restricted number of vpIiipIp^ 
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CHAPTER 5 

CONCLUSIONS AND AVENUES FOR FUTURE WORK 


5.1 Conclusions : 

In this dissertation a vehicle routing and scheduling problem in a multi-product distributed 
manufacturing system is modelled as a standard pickup-delivery problem with time 
windows and additional product sequence constraint. We formulated this problem as a 
mixed integer program and discussed the complexity and inability of the exact methods to 
solve the practical size problems We proposed two approximate algorithms based on 
sequential and parallel insertion procedures and tested few despatch rules for assigning the 
tasks to the vehicles. We performed computational study by designing various types of 
problems. 

The computational study shows that at the best values of slack cost and vehicle waiting 
cost the Parallel Insertion Heuristic is performing always better than Sequential Insertion 
Heuristic. Among the despatch rules applied, the First Come First Serve rule applied to 
Parallel Insertion procedure is performing better in almost all of the cases. In some of the 
cases Systematic Slack Costing rule applied to Sequential Insertion procedure is 
performing well. In terms of average number of vehicles used the First Come First Serve 
despatch rule is performing well for both the algorithms. 
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5.2 Avenues of future work: 


In this work, we assumed there is no capacity restriction for the process centres. One can 
extend this work, by considering capacity restrictions in the process centres. Advanced 
search techniques like genetic algorithms and artificial intelligence techniques can also be 
tried on this problem. 



APPENDIX-A 


PROCESS DATA FACTORS COMBINATIONS WHILE GENERATING 

PROBLEMS 



Process utilisation 

Pr.time/Tr.time 

Slack available 

1 

Even 

High 


2 

Even 

High 

Medium 

3 

Even 

High 

Low 

4 

Even 

Medium 

High 

5 

Even 

Medium 

Medium 

6 

Even 

Medium 

Low 

7 

Even 

Low 

High 

8 

Even 

Low 

Medium 

9 

Even 

Low 

Low 

10 

High-Medium-Low 

High 

High 

m 

High-Medium-Low 

High 

Medium 

12 

High-Medium-Low 

High 

Low 

13 

High-Medium-Low 

Medium 

High 

14 

High-Medium-Low 

Medium 

Medium 

15 

Hi gh-Medium-Low 

Medium 

Low 


16 High-Medium-Low Low High 


17 

High-Medium-Low 

Low 

Medium 

18 

High-Medium-Low 

Low 

Low 

19 

High-Low 

High 

High 

20 

High-Low 

High 

Medium 


2 1 High-Low High Low 


22 

High-Low 

Medium 

High 

23 

High-Low 

Medium 

Medium 


24 


High-Low 


Medium 


Low 


25 


High-Low 


Low 


High 


26 


High-Low 


Low 


Medium 


27 


High-Low 


Low 


Low 














































































APPENDIX - B 


PROBLEM SIZE FACTORS COMBINATION WHILE 
GENERATING PROBLEMS 




No. of operations per product 

No. of process centres 

n 

High 

High 

High 

B 

High 

High 

Low 

3 

High _ 

Low 

High 

4 

High 

Low 

Low 

5 

Low 

High 

High 

6 

Low 

High 

Low 

7 

Low 

Low 

High 

8 

Low 

Low 

Low 
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